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We propose simple Monte Carlo method for calculating parton distri- 
bution in nuclei. Only events satisfying the exact kinematical constrains 
of the corresponding deep-inelastic reaction probing given nuclear distri- 
bution are selected to form the final distribution we are looking for. The 
EMC effect is automatically included by means of two parameters, which 
characterize the change of the nuclear pion field. Good agreement with 
experimental data in the broad range of variable x is obtained. 

PACS numbers: 12.38Aw; 12.38.Lg; 12.39.-x 

1. Introduction 

The study of partonic distributions inside the nucleon and nuclei has 
already long history [1]. Here we shall describe simple Monte Carlo method 
for calculating parton distributions in nuclei with special emphasis on the 
observed differences between such distributions for free and nuclear nucleons 
(known under the name EMC-effect) [2]. The parton picture of nucleon was 
originally formulated in the infinite momentum frame [3]. However, nuclear 
effects are more visible in the nucleon (or nuclear) target rest frame. Such 
change of frames has profound dynamical consequences. Whereas in the 
infinite momentum frame partons can be treated as on-shell objects (with 
some small current masses), in the rest frame they are dressed, far off-shell 
objects with masses consisting substantial part of the nucleon mass. The 
important point is how to find in this frame the proper parton energy. The 
nucleon is no longer contracted by Lorentz transformation and the corre- 
sponding interaction picture is complicated one. The nuclear parton density 
distributions are therefore obtained by generating initial parton momenta 
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in nucleons and calculating the corresponding light-cone longitudinal mo- 
mentum fractions, the so called Bjorken x = xb, which must be identical in 
both frames. 

Our approach is based on the model where valence parton momenta in 
hadron at rest are calculated from a spherically symmetric Gaussian dis- 
tribution with a width derived from the Heisenberg uncertainty relation, 
whereas the sea parton contributions result from similar gaussian distribu- 
tion but with a width dictated by the presence of virtual pions in hadron [2] . 
When going to the nuclear case these initial gaussian momentum distribu- 
tions are changed accordingly in order to account for the presence of nuclear 
medium (like rescattering effects or changes in the virtual pion clouds in the 
nuclear matter). The energy momentum conservation is always strictly im- 
posed and plays vital role in getting our results. The nuclear parton density 
distributions are then obtained by generating initial parton momenta of nu- 
cleons and calculating the corresponding light-cone longitudinal momentum 
fractions Xb- 



2. Deep Inelastic Scattering 

We start with short recollection of necessary theoretical points. In the 
deep inelastic electron-nucleon scattering the cross section is given by: 

da = L„ U W^ (1) 

where the hadronic tensor W^ v can be expressed in terms of electromagnetic 
currents as: 

WW = J Aexp(^e) < P | <V,(0) | p > . (2) 

The spatial variable £ is directly connected to the correlation length existing 
in this process. The virtual foton momentum transfer is given by: 

q=(v,0,0,-y/" 2 + Q 2 ). (3) 

In the Bjorken limit Q 2 — > oo the x = Q 2 /2Mu ~ j + /p + is fixed and 
q 2 jv 2 — > 0. In this limit q~ = q° — q 3 — > oo but q + = —Mx/y/2 remains 
finite. These imply £ + — ► and £~ < y/2/Mx. Alltogether this gives the 
following restrictions for £: 

£o < 1/Mx, t z = z< l/Mx. (4) 

We have therefore two resolutions scales in deep inelastic scattering: 

(i) l/\/Q 2 connected with virtuality of 7 probe. Any two different Q 2 
resolutions are connected via well known A-P evolution equation [4] . 
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(ii) = z being distance how far can propagate the anti-quark in the 
medium, see Fig.l. Notice that small x means a relatively large cor- 
relation length z. Because final state quark interaction within the 
nuclear envirovment is practically not known the small x region opens 
room for different phenomenological models and in this paper we shall 
propose a new mechanism for the nuclear shadowing. 

Deep inelastic scattering of electrons on nuclear targets can be regarded 
as two step process: at first nucleus is replaced by composition of nucleons 
and (effective) pions representing quanta of nuclear binding forces, then im- 
pinging electrons interact with partons (quarks) composing those nucleons 
and pions. Formally it means that nuclear partonic distribution (structure 
function) F A can be written as convolution, 

F A (x A )/x A = Aj J dy A dx5(x A -y A x)p A (y A )F 2 N (x)/x, (5) 

of nucleon distribution function in the nucleus, p A {y A ), and structure func- 
tion of free nucleon, F 2 (x). The x A /A is the ratio of quark and the nu- 
cleus longitudinal momenta (i.e., it is the Bjorken variable for the nucleus) 
whereas y A /A = p + jP\ is the ratio of the nucleon and nuclear longitudi- 
nal momenta. Finally x is the ratio of the quark and nucleon longitudinal 
momenta (the Bjorken variable for the nucleon, x = xb)- 

In the on-shell relativistic approach the nucleon distribution function is 
connected to the well known nucleon spectral function Sn (discussed in the 
next section): 

p A (y A ) = J d A p5 (y A /A - p+/P + ) S N (p°, p). (6) 

In the mean field approximation Sn(p°,p) = n(p)5{p — {m + e(p)), where 
e(p) is the nucleon single particle energy. This expression should be cor- 
rected by the incident flux factor (see eq.(8). It turns out, however, that 
even then one cannot describe the EMC data without inclusion of higher 
order nucleon-nucleon correlations (with additional free parameters and all 
uncertainties of off shell behavior of nucleons it brings in) [5]. 

3. Nuclear Relativistic Mean Field. 

In the nuclear relativistic mean field (RMF) method electrons collide 
with nucleons which are moving in some constant average scalar and vector 
potentials in the rest frame of the nucleus according to the equation: 



[a ■ p + P(m + U s ) - (e N - U v )]ip = 0. 



(7) 
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Here U s = -g 2 /m 2 p s , U v = V^S^o = 9v/ m lp with 9i, m i (* = s ,v) being 
the scalar or vector meson coupling constants and their masses, respectively, 
whereas p s = J2i tptfi^i an< ^ P = Si ^t^i are ^ ne sca l ar an d the fourth com- 
ponent vector densities, respectively. The scalar and vector mean fields were 
investigated successfully in nuclear Dirac phenomenology with the values 
U v = 300, U s = -400 [MeV(p/p )] Po = 0.17 fm^ [6]. It turns out that 
Us and Uy usually cancel each other in the energy or external response 
functions but their relatively big values can explain the enhancement of the 
spin-orbit part of nucleon-nucleus optical potential 1 . It was shown [8] that 
P A (va) depends on both the scalar and vector nuclear fields: 

P A (VA) = -J j^S N (p°,p)(l+p* 3 /E*(p))6(y - (p„ +pS)/ M ), (8) 

where p is equal to the nucleon chemical potential, factor (1 + p\/E (p)) 
corrects (6) for relativistic effects and the nucleon spectral function is taken 
in the RMF approach to be equal S N = n(p)5(p° - (E* (p) + U v )). Eq. (8) 
can be simplified and written as: 

p A {y A ) = \{v\-{y A -l) 2 )/v\, (9) 

where vf = (p/m)(pp / 'E* F ), v\ = Pf/E* f and y is restricted to region 
< (E F — pf) < my < (E F +Pf)- It means that all nuclear dependence 
is hidden in the nucleon chemical potential, which is, however, too weak 
(about 8 MeV smaller than m) to reproduce the minimu seen in the EMC 
data [9] at x = 0.7 (cf. curve (a) in Fig. 2). 

4. Proposed model for parton distribution in nuclei 

Let j denote four-momenta of struck parton (probed by current with 
virtuality Qq) selected (for valence quarks) from Gaussian primordial dis- 
tribution with width 0.172 GeV, and r the respective four-momentum of 
hadronic remnants. Let also W and W denote their respective invariant 
masses. Events are accepted if: 

< j 2 < W 2 , < r 2 < W' 2 (10) 

The sea parton distribution is given by the convolution of the pionic compo- 
nent of the nucleon, f 1T {x; Qq), and the parton structure of pion, f p i on {x; Qo), 

1 One of the advantages of RMF approach is the equation of state it leads to. For 
example, in Walecka model [7] it gives properly the saturation point for nuclear 
matter and no density saturation for neutron matter. 
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obtained from the same Gaussian primodial distribution as used for valence 
partons 2 , 



f s (x; Qq) 



~~fir(yi Qo)fpion(x/y, Qq). 



(11) 
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Fig. 1. 



Fig. 2. 



The characteristic behaviour of the sea partons is derived from the pion 
distribution in the nucleon, which was again parametrized by Gaussian dis- 
tribution, this time with a smaller width equal to 0.052 GeV. The overall 
partonic distribution is therefore given by (Qq = 1 GeV 2 ) 



F 2 (x;Q 2 ) 



Mx,Ql) + f 8 (x;Ql) 



(12) 



Our results for H(x)=F2 e / are presented with Monte Carlo error as 
(b) in Fig. 2. For small x the crucial factor turns out to be the change of 
nuclear virtual pion cloud connected with exchanged mesons responsible for 
the the nuclear forces [2]. In order to be able to fit data in this region we 
have to adjust in our model the value of the parameter which determines the 
relative number of the (effective) intermediate pions (assumed to mediate 
nucleon- nucleon interaction). In our model only part of pions contribute 
to the sea quark structure function of the nucleon whereas the other part 
is responsible for the the nucleon-nucleon interaction. Because for small x 



Actually, in this way we obtain the light cone target rest frame variable x = xlc = 
k + /p + (where p = (M, 0, 0, 0)) with a fixed resolution Q 2 — Qq, whereas experi- 
mentally accessible is Bjorken variable x = XBj = Q 2 /(2p ■ q). However, in the ratio 
presented in Fig. 2, the corrections introduced when processing from one x to another 
cancell and therefore the experimental x in Fig. 2 will be identified with our x. 
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the scale z shown in Fig. 1 becomes comparable or bigger then the nucleon 
size, the expected range of nuclear forces also grows accordingly. Therefore 
for small values of x the part of pionic contribution can "disappear" during 
the interaction with the electromagnetic probe and consequently gives no 
contribution to the nuclear structure function. In our case up to 12% of 
pions are excluded from interaction by this mechanism. 

For intermediate x (~ 0.7) the minimum in the EMC ratio is obtain 
by adjusting the nucleon size in the medium (strictly speaking the size of 
the valence parton momentum distribution). The corresponding increase 
of this width from 0.18GeV to .172GeV produces both the minimum for 
~ 0.7 and the maximum for x around 0.1. In hadronic language this change 
corresponds to some spreading of the pionic cloud outside the nucleon 3 . 

5. Summary 

We obtain very good fit to the data on deep inelastic scatterings with 
56 -Fe using only two physically motivated parameters 4 . The first describes 
decrease of width of the primordial gaussian valence quark distributions 
(it points towards the possible deconfinement of quarks in nuclear matter 
and to chiral symmetry restoration). The second parameter diminish the 
amount of nuclear pions below x = 0.1 due to the shadowing effect. 
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